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Abstract. One-dimensional unitary scattering controlled by non-Hermitian (typically, 
■pT-symmetric) quantum Hamiltonians H ^ H* is considered. Treating these operators 
via Runge-Kutta approximation, our three-Hilbcrt-spacc formulation of quantum theory 
is reviewed as explaining the unitarity of scattering. Our recent paper on bound states 



Znojil M., SIGMA 5 (2009), 001, 19 pages[ [arXiv:0901.0700| is complemented by the text 



on scattering. An elementary example illustrates the feasibility of the resulting innovative 
theoretical recipe. A new family of the so called quasilocal inner products in Hilbert space 
is found to exist. Constructively, these products are all described in terms of certain non- 
equivalent short-range metric operators 9^J represented, in Runge-Kutta approximation, 
by (2R — l)-diagonal matrices. 

Key words: cryptohermitian observables; unitary scattering; Runge-Kutta discretization; 
quasilocal metric operators 
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1 Introduction and summary 

In paper I p] we summarized the situation in which certain very complicated quantum bound- 
state Hamiltonians f) = t)( p ) = fjt (acting in the usual physical Hilbert space were as- 
sumed solvable via a Dyson-inspired transition to their non-Hermitian but perceivably simpler 
isospectral partners H = S7 -1 rjf2 ^ W. The latter auxiliary operators were assumed acting in 
a doublet of the first auxiliary and the second auxiliary Hilbert spaces H^ F ' S \ respectively. The 
latter space 

W (S) 

was finally assumed endowed with an unusual inner product (a,b)^ defined 
by the formula (a, b)^ := (a, 06)( F ) in terms of the usual Dirac's inner product (a, b)^ and of 
the so called metric operator O = 0t = f^il > 0. 

In the present continuation of paper I we shall extend our three-Hilbert-space formulation of 
quantum theory to a class of models of scattering. For illustrative purposes we shall discretize 
the axis of coordinates into Runge-Kutta lattice [2]. This discretization is not a mandatory 
ingredient of the approach but its use will facilitate our explicit constructions. In particular, 
it will enable us to demonstrate that for a selected sample Hamiltonian H one can construct 
several alternative metrics which are all compatible with the unitarity of the scattering-state 
solutions of the underlying Schrodiner equation. This quantification of the ambiguity of will 
be a core of our message to physicists. Unexpectedly, our unitarity-supporting illustrative 0s 
will prove "quasilocal", i.e., they will emerge as short-ranged continuous-coordinate limits of 
(2R — l)-diagonal Runge-Kutta approximation matrices = Qr with integer R = R(h) such 

that lim hR(h) = 0. This will be our main new mathematical result. 

h->0 



*This paper is a contribution to the Proceedings of the 5-th Microconference "Analytic and Algebraic Me- 
thods V" . The full collection is available at |http: / / www.emis.de /journals/SIGMA /Prague2009.html | 
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The detailed presentation of this result will be preceded by a short review of motivation 
(Section [2D and of the currently well-established theoretical understanding of non-Hermitian 
models of bound states (Section [3|). The current status of extension of this version of quantum 
theory to scattering scenario will be summarized in Section HJ 

In Section [5] we pick up one of the elementary concrete models [3] and use it to outline our 
method applicable in the dynamical regime of scattering. The core of our message is formulated 
in this section. We fix there the level of approximation given by the Runge-Kutta lattice-spacing 
h > and describe the construction of metric operators O = Q(H). For our toy model H, in 
particular, we list the first few explicit samples of B = Qr(H) possessing the {2R — l)-diagonal 
matrix structure. These matrices are constructed by computerized symbolic manipulations at 
R = 1,2, ... ,7 '. Our knowledge of these solutions enables us to conjecture (and, subsequently, 
prove) formula which defines all the eligible quasilocal metrics Or(H) at any R in closed form. 

In complementary and concluding Sections [6] and [7J we construct underlying Dyson opera- 
tors f2 and discuss some of their potentially most relevant descriptive properties as well as several 
possible consequences of their use in model-building. 



2 Context: T'^-symmetric Hamiltonians 

Among all of the simplified phenomenological models used in quantum mechanics a prominent 
role is played by the families described by Schrodinger equations of ordinary differential form 
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■9(x) + V(x)V(x) = EV(x). (1) 



Centrifugal coefficient L(L + 1) may be admitted to vanish. Variable x is then interpreted as the 
one-dimensional coordinate. Requirement \&(x) G L2(R) is being imposed upon bound states 
while the scattering solutions are characterized by asymmetric asymptotic boundary conditions 
at all real k = \J E — V(oo), 

( Akx i TJp—iKX rr ^ _1 

*(x) = \t iK Z ' x<< /' (2) 

v ) y Te lKX , x > 1. w 

At nonvanishing centrifugal coefficients one demands that L{L + 1) > —1/4 while variable x 
is treated as the radial coordinate in D dimensions, x = \r\ € (0, oo). An additional, specific 
boundary condition must be then imposed in the origin in order to guarantee the regularity of 
bound-state as well as scattering solutions Of course, physics represented by equation (P) 
on half line depends again on the asymptotic behavior of potential V(x). 

Usually one distinguishes between the confining, bound-state regime (where V(x) is chosen 
smooth, real and very large at large x — > ±oo) and the scattering scenario (where, typically, 
V(±oo) = and where wave functions ty(x) are selected as free waves at large x — > ±oo). 
Recently, several authors emphasized the existence of another, unusual quantization of the or- 
dinary differential Schrodinger equation ([1]) which is assumed integrated along a non-standard, 
complex contour of "pseudocoordinates" x = x(s) £ C where s S R is a parameter (cf., e.g., 
Sibuya [5] or Bender and Turbiner [6]). 

The early attempts in this direction were treated and accepted as a mere mathematical 
curiosity. For illustration we may recall paper [7] by Buslaev and Grecchi who studied anomalous 
potentials V(x) ~ — x 4 + 0(x 3 ). Equation ([T]) with these "wrong-sign" potentials has been 
considered along a shifted line of x = x^ BG ^ (s) = s — ie. By this trick, the potential (which looks 
asymptotically repulsive) has been made manifestly confining. In parallel, the usual requirement 
of Hermiticity has been replaced by 'PT-symmetry, i.e., by left-right symmetry in the complex 
plane of x (cf. also [8] for more comments). 
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During the last ten years we witnessed a quick growth of interest in 'PT-symmetric models. 
Within physics community this interest has been inspired by Bender's and Boettcher's influential 
letter [9] where the authors emphasized the deeply physical appeal of the combination of the 
parity-reversal symmetry (mediated by the operator V) with the time-reversal symmetry (cf. 
operator T). They conjectured and demonstrated, by approximate methods, the reality (i.e., 
in principle, measurability) of the bound-state spectra for the whole family of specific VT- 
symmetric toy potentials 



A few years later the validity of the conjecture has rigorously been proved by Dorey, Dunning and 
Tateo |10] and by Shin [11]. This clarified, for bound states, the mathematics of complexification 
of the argument x of wave functions ip( x ) and/or of potential V(x). 

The questions emerging in scattering regime have temporarily been left open. 

3 New horizons: cryptohermitian observables 

The success of the latter generalization of the class of quantum bound-state models opened 
several new theoretical directions of research since in many concrete examples the complex 
value of "coordinate" x ceased to be tractable as a measured or measurable quantity. For this 
reason the theoretical as well as practical acceptance of the internal mathematical consistency of 
apparently non-Hermitian models exemplified by equations (H|) + © took some time in physics 
community (cf., e.g., review paper [12] for historical comments). The paradoxes seem to be 
clarified at present and, for bound states, a new pattern is established for implementation of 
the formalism of textbook quantum mechanics. 

The key to the safe return from equations (fl|) + (and the like) to textbooks can be seen 
in the concept of cryptohermiticity (the word invented, very recently, by Smilga [13]) which, in 
essence, means that an operator H which appears non-Hermitian in Hilbert space 7Y (F) (where 
the superscript indicates the (user-)friendliness of the most current mathematical definition of 
the inner product - see below) may be reinterpreted, under certain circumstances, as safely 
Hermitian in another, Hilbert space (where the superscript stands for "standard" physics). 
The third Hilbert space emerges, quite naturally, as a space which is unitarily equivalent 
to TiS^ (i.e., it represents strictly the same physics so that for the purposes of physical predictions 
we have an entirely free choice between and TC^) but which is the only space encountered 
in textbooks (i.e., the metric operator in TL^ remains conventional, identically equal to the 
unit operator). 

One of the first papers presenting and reviewing such an abstract idea in a concrete applica- 
tion in nuclear physics appeared more than fifteen years ago [13] • The study of heavier atomic 
nuclei has been shown simplified by the mapping of the usual, complicated physical Hilbert space 
of states 7~^ p ) (describing fermions and possessing, therefore, complicated antisymmetrization 
features) upon another, bosonic (i.e., manifestly unphysical) auxiliary space TL^\ This made 
the calculations perceivably facilitated. Whenever needed, the return to fermionic wave func- 
tions mediated by the not too complicated Dyson mapping 0, remained feasible via a unitary 
equivalence between and (more details may be found in pQ). 

The combination of a comparatively narrow applicability in nuclear physics with a rather 
unusual mathematics caused that the methodical appeal of the non-unitary-mapping idea re- 
mained virtually unnoticed until its re-emergence in 'PT-symmetric context. At present we wit- 
ness a massive revival of interest in the parallel use of the three alternative representations TL^\ 

and of the same physical quantum system. The invertible Dyson-type mapping 

relates the spaces and Ti.^ as well as the respective Hamiltonians or operators of other 




5 > 0. 



(3) 



observables. 
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The manifest violation of the unitarity by the mapping ft ^ (f^) -1 requires due care in 
applications. For example, the physical Hamiltonian f) must be Hermitian (i.e., more strictly, 
essentially self-adjoint in the physical Hilbert space TC^ exemplified by the fermionic space in 
the above-mentioned concrete application). It may also be perceived as a transform of another 
operator defined as acting in another vector space, f) = QHQ" 1 . The latter "Hamiltonian" F 
is, by construction, defined and manifestly non-Hermitian in space Tt^ where, trivially, 



II 



-i\t 



1\ T 



OF©' 



(4) 



We abbreviated © = calling this new operator, conveniently, a metric. It is defined as 
acting on the kets in both Hilbert spaces TC^ and TC^ . This mathematical ambiguity proves 
inessential in applications and it is also easily clarified by the following elementary graphical 
pattern using notation of paper I, 



map f2 / 
math. OK in 1W> 
- computable 



physics OK in 
ket | tp y = uncomputable 

~7 \~~ 



map Qfl 1 =I 



map SI 1 
all OK in 

/A = the same 



This diagram reconfirms that the Hermiticity of the physical Hamiltonian f) = QHQ, 1 = f)t = 
acting in the physical Hilbert space is equivalent to the cryptohermiticity 

constraint ^ imposed upon its partner H = / Ft in H^l 



3.1 The ambiguity of the metric © = @(H) 

Hesitations may emerge when one imagines that the Hermitian conjugation itself is not a unique 
operation and that it may be altered [15]. In this way equation ([4]) may be interpreted as assign- 
ing several alternative physical metric operators to single Hamiltonian, O € (©i(F), 02(F), . . .). 
The first thorough discussion and clarification of this problem of ambiguity has been published 
by Scholtz et al. [14]. They emphasized that besides the Hamiltonian F itself, any other op- 
erator O = 0\,0 = O2, ■ ■ ■ of an observable quantity must obey the same cryptohermiticity 
relation as H = Oo itself. In opposite direction, for a given set of observables in 7i^\ each 
eligible metric operator O must remain compatible with all of them, 

QO j = 0}e, j = (0,)l,2,...,j max . (5) 

These requirements reduce the ambiguity of O so that at a suitable integer j max the metric 
O = 0(F) may become, in principle, unique. 

The ambiguity of O may acquire an enormous theoretical importance as well as phenomeno- 
logical relevance. The abstract understanding of the ambiguity of Dyson maps did already 
prove crucial not only for the above-mentioned efficient description of bound states in nuclear 
physics but also in a clarification of the role of the charge or ghosts in field theory [16], etc. 



3.2 Classical limit in the case of cryptohermitian Hamiltonians 

Persuasive demonstration of relevance of the apparently purely mathematical subtlety of the am- 
biguity of metric can be mediated by the popular Bessis' and Zinn- Justin's (BZJ, [T7]) imaginary 
cubic oscillator 

H^P) = E n \^), F( BZJ ) = --^ + i^ 3 , n = 0,l,... (6) 
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for which the quantum bound-state problem is formulated and solved in Hilbert space 

W (F) = 

L 2 (M) and for which the spectrum is real, discrete and bounded below despite the BZJ Hamil- 
tonian being manifestly non-Hermitian. 

Routinely, the Dyson mapping may be employed to convert the elements \ip^) £ H.^ of 
the above-mentioned space into their images £ Ti.^ which lie in the correct representa- 

tion Ti,^ of the abstract, textbook Hilbert space of states. Superscript ( p ) stands for "physical" 
and we have tt : -> H (p) and |Vj (p) ) = ST2|V> (F) )- It is the isospectral partner f)( BZJ ) = 

our original Hamiltonian which becomes, by construction, safely Hermitian in 
physical space . 

There is a price to be paid for the clarification of theoretical concepts. In our illustrative 
example ([6]) it was difficult to prove that the spectrum of the bound-state energies is all real |10j . 
In [18] one finds that even in the weak-coupling regime with a very small value of g = e one 
arrives at an impressively complicated physical representative of Hamiltonian in 7i^ p \ 



^(b ZJ ) = ntfCBz^n-i = t + *_ (| x6) J_| + 22 J_ j + (510 + 10V J j, i 

+" — V'"" 1 " Ui I* 3 . e 



p s 3 x [ ' p 5 } ^) 



p 6 



(7) 



This formula contains parity operator V and anticommutators {•,•}, it uses the Fourier-trans- 
formed multiplication-operator representation of the powers of p in L 2 (R) and, finally, it varies 
freely with four real parameters \j and Kj, j = 1, 2. Persuasively, formula ([7} demonstrates that 
the explicit form of true Hamiltonians f) = f)^ may become threateningly or even prohibitively 
complicated. One cannot be expected to perform practical calculations using this representation 
of the Hamiltonian. 

The latter remark seems to imply that the explicit knowledge of maps H <-> f) or, for other 
observables, A <-> a may be skipped as redundant. Actually, this is just a partial truth since 
operators f) or 0, may remain unavoidable, e.g., during the analysis of time-dependent systems 
(cf. [19] and also the recent preprint [20J in this respect). Moreover, in model ([7]) one could 
select A2 = k-2 = in order to preserve 'PT-symmetry. Last but not least, the "unfriendly" 
Hilbert space opens the way towards the classical limit of the system. From the imaginary 
cubic example ([7]) one arrives at the classical Hamiltonian function h c defined in the classical 
single-particle phase space of the coordinate p c and momentum p c , 

h^ = ^ + e* 9 (p c )4 + 0(e% *<Pc) = g (8) 

(cf. equation (62) in [H]). Thus, the limiting transition h — > reveals a hidden sextic-oscillator 
nature of the imaginary cubic forces in the weak-coupling regime. 



3.3 Extended Dirac's bra/ket notation 

The popular emphasis on the non-Hermiticity of simplified, artificial H 7^ in should 
be perceived as attracting attention but slightly misleading. It reflects just the mathematical 
property (implicitly, of the only relevant, correct and physical operator fj) with no direct con- 
nection, say, to the principle of correspondence (cf. the preceding paragraph). In contrast to 
the descriptions of open systems [21] where the probability is not conserved the calculations 
of probabilities must all be performed inside Ti^ or in one of its unitarily equivalent partner 
Hilbert spaces Ti^ where 6//. We strongly recommend an explicit reference to the space in 
which one resides, because 



6 



M. Znojil 



• the apparent and misleading conflict between the absence of the correct physics in TC^ 
and the lack of the computational feasibility in is naturally resolved by the metric- 
dependent description of the system in the standardized Hilbert space 

• the Dirac notation can be used after a graphical adaptation of the bra- and ket-symbols 
when one speaks about the elements of the formally different Hilbert spaces 
andft( p ). 

More explicitly, the same state of a quantum system will be characterized by a spiked ket 
symbol in P-space, \ip>~€z 7i^ p \ and by the standard ket symbol in the other two Hilbert spaces, 
|?/>) 6 W( F ' S ). The two spaces are chosen as identical when only their ket elements are 

being considered. These two spaces are identical as vector spaces which should be denoted by 
a dedicated symbol = V^ s ^. They only differ in the respective definition of their respective 
linear functionals, i.e., of the bra-vector elements of the dual vector spaces marked by the prime, 
y( F ) y( s ) . The correspondence between the three dual vector spaces can be represented by 
the following diagram, 



map 



<M g (y( p >)' 



g (y( F >y 



map e=rtn^i 



map f2 



e (y( s >)' 



In this notation it is easy to define {{^f\ = (\l/|0 or to work with the map 0, : \ifj n ) — > \tfj n y. It 
is equally easy to verify the unitarity of the map between the Hilbert spaces W( p ) and 

^\cj»= ((t/#). 

As long as (V^ F ')' 7^ (V^ s ^)' we must distinguish between the explicit definitions of the respective 
operations of Hermitian conjugation. The operation marked by single-cross superscript ' in TL 1 ^' 
must be distinguished from the one active in H {s) and marked by double-cross *. 

The acceptance of the notation recommended in this section enables us to introduce a modi- 
fied bra symbol in 7i^\ (1^))* = ((^| £ (V (s) )' while keeping the traditional (| x I / )) t = (^| £ 
(y( F ))' unchanged. There is no doubt that the use of three Hilbert spaces 

ft(P,F,S) doeg not 

contradict any principle of textbook quantum mechanics. Indeed, all the physical operators of 
observables remain Hermitian in Ti.^ while the other two spaces are just auxiliary. 



4 Scattering and problems with long-ranged non-Hermiticities 

The variability of 0s and/or Sis could inspire optimistic expectations concerning the extension 
of the cryptohermitian description of physics to scattering scenario. Unfortunately, several se- 
rious obstacles were found and formulated by Jones [22]. He argued that Dyson maps O must 
be necessarily long-ranged leading to an apparent violation of causality in scattering (cf. Sec- 
tion SJ] below) . A return to optimism has been initiated in [23] where some of the most striking 
difficulties were circumvented via a replacement of the current, strictly local interactions V(x) by 
their minimally nonlocal alternatives (cf. Section [4.21 below). The first fully satisfactory model 
of scattering has been found in [3] (cf. Section T4.3I below) . A continuation of these developments 
will be reported here immediately after a brief summary of older results. 

4.1 The first problem: The causality-violating waves in the scattering 

In numerous recent applications of the above-mentioned three-space representation of quantum 
theory to bound states a conflict was encountered between the friendliness of the calculations in 
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unphysical and the unpleasant complications arising during the transition to the two alter- 
native physical spaces 7i^ S ' F \ Many concrete Hamiltonians H happened to possess a particularly 
simple form in space and vice versa. 

The transition to the correct space may be perceived as a fairly unnatural operation. The 
more so if one tries to describe the non-Hermitian scattering 125] . For these reasons, the 
counterintuitive character of the textbook Hermiticity of H in unfriendly is reflected by 
the terminology in which the cryptohermitian operators are nicknamed quasi-Hermitian rather 
than Hermitian [141 [26] . 

The main difference between the two Hilbert spaces concerns their inner products 

between elements and \4>). In the most usual coordinate representation language we have 



(V#) (F) = / ^*{x)(j){x)dx = (V#> 



m 



and 



(S) 



^*(x)@(x, x')(fi(x') dxdx' 



lei 



m 



(9) 



(10) 



where Q(x,x') ^ 5(x — x'). The consequences of the transition from equation ([9]) to equa- 
tion (|10p are nontrivial. For example, according to Mostafazadeh [27] the assumption of having 
a short-range and strictly local V = V(x) leads immediately to the necessity of using a strongly 
nondiagonal metric kernel Q(x,x') in equation (|10p . 

The latter long-range nonlocality is easily visualized in Runge-Kutta picture where the co- 
ordinates are represented by a lattice of discrete points x = Xj, j = . . . , —1, 0, 1, . . . and where 
all the differences Xj+i — Xj = h are the same. In this approximation the typical doubly infinite 
matrix @jj = @(xi,Xj) is dominated by the unperturbed unit matrix 0^ = 5ij. The non- 
negligible correction given, e.g., by equation (16) in [22] will be strongly non-diagonal. The 

metric can comfortably be expanded in the sum of elementary matrices, Ylk e ~^ hk ®ij where 
each coefficient is sparse and strongly non-diagonal, 



Q (i) 



1 

1 

1 


1 

1 

1 


1 

1 

1 


1 

1 

1 



1 

1 1 
1 


1 

1 1 
1 


1 

1 1 
1 


1 

1 1 

1 



, .... (11) 



These cross-shaped matrices couple remote coordinates Xi and Xj. A conflict emerges between 
our use of the common Hilbert space TL^ (which keeps trace of the intuitive principle of cor- 
respondence) and of its amendment 7i^ s ' (only there the correct probabilistic interpretation of 
the observables is achieved). The core of the problem lies in the manifest loss of the concept of 
the asymptotically free motion. In words of [25] one has to conclude that the causality- violating 
nature of metric is certainly "changing the physical picture drastically". 



4.2 The second problem: Spatial symmetry violations 

The first steps towards the desirable return to the standard picture of scattering induced by 
the Hamiltonian H which proves non-Hermitian in space were proposed in our com- 

ment [23]. We tried to resolve there those of the paradoxes covered by paper [25] which 
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were available prior to publication [28]. In particular we paid attention to the apparently 
unavoidable presence of causality-violating waves which seemed to emerge in the right spa- 
tial infinity. In the context of the specific non-Hermitian delta-function scattering models this 
violation of causality seems to result from the strict locality of V = V(x). We concluded 
that such an assumption proves too strong and that it must be weakened for the given pur- 
pose. 

In a detailed discussion of the causality- violation paradox we employed the same discretiza- 
tion. We assumed that in Hamiltonians H = —d 2 /dx 2 + V the potential of any form must be 
combined with the tridiagonal Runge-Kutta version of the kinetic-energy operator, 



H = — A + V, 



-A 



2 -1 
-1 2 -1 

-1 2 -1 

-1 2 



(12) 



This helped us to simplify some technicalities and to clarify the way in which the asymptotically 
non-vanishing nondiagonality in 0jj is generated by the diagonality of potentials V = V{xi). 
We then decided to consider only such families of Hamiltonians for which the unpleasant non- 
diagonality of the metric disappears, 



@(xj,yk) « c(xj)Sj ;k , 



>1, 



(13) 



(cf. equation (23) in [23J). Then we only had to recollect that the (crypto) Hermiticity of H 
m is a condition which acquires the utterly elementary form @ in . As a linear- 
equation constraint imposed upon the matrix elements of metric O this equation can serve 
as an independent, nonperturbative source of information about the nonlocalities in met- 
rics ([IT]). 

In [23] we decided to shorten the range of the influence of the non-Hermiticity. After the 
insertion of ansatz (|13p in equation @ a series of our algebraic trial and error experiments 
revealed that in order to achieve a certain internal consistency of our requirements we might 
replace the usual complex and diagonal non-Hermitian matrix V(x{) in H = — A + V(x) by its 
two-diagonal real and non-Hermitian analogue of the form 



y(a,b,c,...) 



—c 

c -b 
b 


—a 


a 


-b 

b -c 
c 



(14) 



A multiparametric Schrodinger Hamiltonian H = —A + V^ a ' b ' c '"^ has been found which is 
Hermitian in the Hilbert space Ti^ where the exact metric operator has the following nontrivial 
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but still fully diagonal form, 



0-5 

0-3 

0-i 






Oi 

03 

05 



(15) 



Matrix elements are explicitly known, 

6 ±1 = (l±a)(l-b 2 )(l-c 2 )(l-d 2 )--- , 
9 ±3 = (l±a)(l±b) 2 (l-c 2 )(l-d 2 )--- , 
6 ±5 = (1 ± o)(l ± bf{l ± c) 2 (l - d 2 ) ■ ■ ■ . 

We separated the "in" and "out" solutions not only in TC^ but also in TC^ . A causality- 
observing physical picture of scattering was demonstrated to exist. Unfortunately, the amend- 
ment of the potentials proved incomplete. We did not manage to remove another shortcoming 
of equation ([15]) where, for a generic set of parameters a, 6, . . ., the asymptotic measure c(xj) — 1 
of the anomaly of the flux in equation (|13p remains long-ranged (cf. a footnote in |25j). 



4.3 Problems resolved: PX-symmetric models of scattering of [3] 



Having accepted the above-cited words of critique we sought for a removal of the long-range 
anomalies. In letter [3] we made another step towards a fundamental theory of scattering based 
on the Hamiltonians which appear non-Hermitian in space TC^ F \ Another, amended family 
of Runge-Kutta Hamiltonians has been found there preserving the kinetic + potential-energy 
structure, H = T + V. We restricted our attention to the following two-diagonal matrix inter- 
action proportional to real coupling parameter g and mimicking the existence of two interaction 
centers at a distance ~ 2M, 

















-9 

9 




























9 

-9 









2M - 3 
columns 








(16) 



The key point was that we achieved an asymptotic spatial symmetry of the metric operator, 
c(xj) = c(-Xj) = 1, \xj\ » 1. (17) 
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The overall picture offered by equation (|17p looks satisfactory, in the discrete-coordinate ap- 
proximation at least. The metric assigned to our Hamiltonian H remains diagonal, even though 
it is still different from the unit matrix. We can reasonably expect that for the majority of 
interactions the diagonality of this matrix may easily happen to lead just to a smooth change 
of the measure in equation © giving just the replacement dx — > dfi(x) in the limit h — > 0. In 
spite of having a nontrivial Hilbert space, the changes due to O ^ I may remain inessential and 
the limiting metric O may be called "local". 

In contrast, whenever we relax the constraint of diagonality and admit a {2R — l)-diagonal 
metric with R > 2, the decrease of the discretization length h — > may lead, in general, 
to a matrix which can be, say, a polynomial function of the kinetic energy operator A of 
equation (|12p . The h — > limit of the discrete inner product would not lead to the mere 
double-integration formula (llOj) but rather to a very complicated (e.g., integro-differential) 
recipe. 

In our present text we intend to omit the discussion of all the similar mathematical subtleties, 
emphasizing only that they might play an important role even in the present class of illustrative 
examples. This is also the reason why we call our band-matrix metrics quasilocal. In the limit 
h — > 0, their range will certainly shrink to zero but the same feature characterizes also the 
kinetic-energy operator A. 

5 Band-matrix metrics: their meaning and construction 

The fact that we can only reconstruct the reality from our discretized scattering models via the 
limiting transition h — > has an important constructive aspect. Let us recollect that the core of 
the consistence of the theory at h = lies in a reconciliation of the non-Hermiticity of differential 
equation (fTJ) in with the asymptotic locality of boundary conditions ([2]). We saw that at 
h > both these problems found their resolution in [3] where H was simply re-interpreted as 
Hermitian in TC^ and where the diagonality of the metric at h > offered an easy guarantee 
of its locality in the limit h = 0. 

Now we have to add that the classification of possibilities as presented in [3] is incomplete 
because we required there the strict locality of the metric. On this background we may briefly 
characterize our present new results as a broadening of the picture of unitary scattering as 
offered in [3]. 

Our present key idea is that the standard form of the asymptotic boundary conditions for 
scattering (defined, naturally, in terms of plane waves) would not be lost even if we assume 
that our metric becomes weakly non-diagonal. For example, certain constant tridiagonal 
matrices constructed at a given lattice constant h > may appear as proportional to 
the Laplace operator (i.e., to the kinetic-energy operator) in the continuous coordinate limit 

This point of view can be generalized to all of the {2R — l)-diagonal matrices of the metric 
where the integer R remains fixed during the decrease of h — > 0. The resulting limits may 
often prove expressible in terms of the powers of the Laplace operator so that we may call all of 
these h = metrics quasilocal. 

As a consequence, a broad family of new metrics could be obtained in the limit h = 0. 
They would still, by construction, leave the scattering unitary. The resulting picture of physics 
will vary with the choice of their structure at h > 0. In order to be able to understand this 
situation we would need to be able to construct the {2R — l)-diagonal matrices in closed form 
at any h > 0. 

This possibility has been ignored before. Let us now fill this gap via our solvable non- 
Hermitian example. 
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5.1 The choice of the model 



Let us select the class of interaction models (fT6|) of [3] and, for illustration purposes, let us pick 
up the special case where, formally, 2M — 3 = — 1, 



H 



2 -1 

-1 2 -1-g 

-1 + g 2 -1 + g 

-1-g 2 -1 

-1 2 



(18) 



This is a tridiagonal, doubly infinite Hamiltonian matrix which can be complemented 
by a discrete version of asymptotic boundary conditions (|2|) for scattering wave func- 
tions. The details including the explicit formulae for the reflection coefficient R(E) 
and for the transmission coefficient T(E) may be found in [3] showing that \T(E)\ 2 + 
\R(E)\ 2 = 1. This means that the standard probability flow is conserved (i.e., the pro- 
cess is unitary) even when the scattering is controlled by our non-Hermitian Hamilto- 
nian lfT8jh 

The unitarity of the model is related to the existence of a metric which is asymptotically 
diagonal. At any nonvanishing h > we may easily find the fully diagonal metric containing 
just a single anomalous matrix element, 



6 = Bi 



(1 + <?)/(! -<?) 



(19) 



Having worked in discrete coordinate representation we may conclude that our matrix ©i com- 
mutes with the usual diagonal operator of the coordinate. The anomalous matrix element 
z\ = (1 + g)/(l — g) = 1 + 2^/(1 — g) degenerates to 1 in the no-interaction limit (i.e., Hermitian 
limit) g — ► 0. One should notice that the spectral singularity [29] emerges at g = ±1. 



5.2 The existence of band-matrix metrics with 2R — 1 < 9 diagonals 

Our choice of the elementary toy-model Hamiltonian (I18p facilitates the constructive search 
for the alternatives to the diagonal metric (|19p . Indeed, whenever we perceive the class 
of admissible metrics O = 0(H) as composed of all the possible Hermitian, positive and 
invertible matrices which are compatible with the Hamiltonian-dependent quasi-Hermiticity 
constraint (|4|) we may (and shall) take this equation simply as a linear algebraic sys- 
tem of (infinitely many) equations which are to be satisfied by our pair of matrices H 
(which is given in advance) and O (which is to be specified as one of solutions of equa- 
tion (HD). 

There exist several technical obstacles which make such a project not entirely trivial. 
Firstly, we have to employ some "manual" linear algebra in order to reduce the infinite- 
dimensional equation to its finite-dimensional subsystem. By the trial and error tech- 
niques we succeeded in revealing that such a reduction may even preserve many parallels 
with diagonal ©i. In this spirit we decided to search for certain doubly- infinite-matrix 
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solutions Or which exhibit a very specific band-matrix property (Qr)^ m = whenever 
\k — m\ > R. Empirically we revealed that each such sparse, (2R — l)-diagonal mat- 
rix solution @r may, furthermore, be required to contain solely R nonvanishing diago- 
nals. This possibility is supported by the tridiagonal matrix ansatz for Qr at R = 2 gi- 
ving 



e 2 



i 

1 l+g 

l+g l+g 

l+g 1 

1 



A reconfirmation is found at R = 3, with 



1 o 



i 

o 



i 
o 

l-9 2 



l + g 

1-9 2 



i + g 



(l+g)(l-2g 2 

1-9 


1+9 



1-9 2 
l + g 



i-g 2 
o 
l 



1 

1 
1 



and at R = 4, with 



0i 



'•• 1 '•• 1 

10 10a 

'•• 1 •■. b ■•■ 

1 b c 



c 6 1 

'■■ b '■■ 1 '■■ 

a 1 1 

1 '•• 1 '■■ 



where we abbreviated a = 1 + g, 6=1 — g 2 and c = (1 + g)(l — 2g 2 ). Finally, we also solved 
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equation with the similar R = 5 ansatz 



e. 



1 i 
i 

i i 

i 

1 b 

a 

b 



d d 
d d 



and we obtained, again, a = 1 + g, b = 1 — g 2 , c = (1 + g)(l — 2g 2 ), d = (1 — g 2 )(l — 2g 2 ) and 
{l+g)(l-2g 2 f 



~5 



1-5 



5.3 The extrapolation and construction of all the metrics ©^ 

As a result of our computer-assisted exact and systematic reconstruction of the metrics B = 0& 
compatible with Hamiltonian H via quasi-Hermiticity condition @ we arrive at the following 
extrapolation of our k < 5 rigorous symbolic-manipulation results to any integer k > 5. 

• By construction, each metric 0^ possesses the form of a doubly infinite and symmetric 
(2k — l)-diagonal real matrix with k nonvanishing diagonals interlaced by k — 1 zero 
diagonals. 

• Up to their central fc-plets, all the matrix elements of along each non-vanishing diagonal 
are equal to one. 

• If we assume k > 2 and abbreviate a = 1 + g and b = 1 — g 2 the leftmost non-trivial fe-plet 
(a, b,b, . . . ,b,b, a) of non-unit matrix elements of our metric contains two boundary a's 
complemented by the (k — 2)-plet of 6's. 

• After we abbreviate c = (1 + g) (l — 2g 2 ) and d = (1 — g 2 )(l — 2g 2 ) and assume that k > 4, 
we obtain the subsequent fc-plet in the explicit form (b, c,d,d, . . . ,d, d, c, b) filled by k — 4 
d's. 

• With e = (1 + g)(l - 2g 2 ) 2 and / = (1 - g 2 )(l - 2g 2 ) 2 and k > 6 we get the next fc-plet 
(b, d, e, /,/,... , /, /, e, d, b) containing k - 6 fa. Etc. 

The general pattern is obvious: the quadruplets of matrix elements sitting in the four (viz., 
"North", "East", "South" and "West") corners may be ordered in an inwards-running sequence 
of vertices a = A\, c = A2, e = A3, ... (cf. Tabled]). These values are complemented by 

the related multiplets of the wedge-filling elements b = B±, d = B2, f = -B3, At the 

odd subscripts there emerges an anomalous, non-polynomial central matrix element z\~ equal to 
fractions a/(l — g), c/(l — g), e/(l — g), . . . at k = 1, 3, 5, . . ., respectively. 
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Table 1. Matrix elements of the metrics 0#. 



R 


central element 


corner element 


wedge element 


1 


z 1 = A x /(l-g) 


— 


— 


3 


z 3 = A 2 /(l - g) 


a = A x = (1 + g) 


b = B 1 = (l-g 2 ) 


5 


z 5 = A 3 /(l - g) 


c = A 2 = (l + g)(l-2g 2 ) 


d = B 2 = (l-g 2 )(l-2g 2 ) 


7 


z 7 = A 4 /(l - g) 


e = A 3 = (l+g)(l-2g^ 


f = B 3 = (l-g 2 )(l-2g 2 ) 2 


9 


z 9 = A 5 /(l - g) 


a 4 = (i + <?)(i- V) 3 


S 4 = (l-5 2 )(l-2 5 2 ) 3 


2k + 1 


z 2 k+i = A k+1 /(1 - g) 


A k = (l + g)(l-2g 2 ) k - 1 


B k = (l-g 2 )(l-2g 2 ) k - 1 



The existence of the above extrapolation rules enables us to formulate a simplified ansatz for 
the next matrix at R = 6, 



i i 
i i 



i 

i i 
i 



i 



i 

b 

d 
d 



1 

a 

b b 

c 

d d 

e 

e e 
e 

d d 
c 

b b 

a 

1 



d 
d 



b 

c a 
b 

b 1 
1 

1 1 
1 

1 1 



It passes the test after insertion in equation @. The tedious computer-assisted direct solution 
becomes replaced by an elementary verification of the choice of e = (1 + g) (l — 2g 2 ) . 

In a climax of our analysis let us now complement our algorithmic set of extrapolation 
rules by their rigorous proof. For this purpose we may select an odd integer R = 2k + 1 
and reinterpret Table [T] as a recurrent pattern (i.e., the proof proceeds by mathematical in- 
duction). The corner element A k is most easily deduced from the previous line since A k = 
(1 — g)z 2 k-i and only the values of B k and z 2 k+i must be deduced from the quasi-Hermiticity 
condition (j3]) rewritten as a doubly infinite matrix set of linear equations using the input 
Hamiltonian (|18p . 



[H^e 2k+1 -G 2k+1 H} mn = 0, m,n = 0,±l,. 



(20) 



Cryptohermitian Picture of Scattering Using Quasilocal Metric Operators 



15 



This task is facilitated by the tridiagonality of H and it may be also guided by our last illustrative 
sparse-matrix ansatz for Qr at the next odd R = 7, 



e 7 = 









1 




1 




















1 




1 




a 
















1 




1 




b 




b 












1 




1 




b 




c 




b 












1 




b 




d 




d 




b 








1 




b 




d 




e 




d 




b 








b 




d 




f 




f 




d 




b 




a 




c 




e 




Z7 




e 




c 




a 




b 




d 




f 


f 




d 




b 








b 




d 




e 




d 




b 




1 








b 




d 




d 




b 




1 












b 




c 




b 




1 




1 












b 




b 




1 




1 
















a 




1 




1 




















1 




1 









Having the latter structure and illustration in mind we can contemplate any R = 2fc + 1 and 
verify, by direct calculation, that the rightmost unknown in the last line of Table Q] (i.e., the 
value of the wedge element Bk) emerges from equation (120p as defined, recurrently, at any 
one of the eight pairs of the subscripts (m, n) = (— 2,±1), (m,n) = (— 1,±2), (m, n) = 
(1,±2) or (m,n) = (2, ±1). Similarly, the leftmost and last unknown of Table Q] (i.e., the 
value of the central matrix element Z2fc+i) becomes determined, in terms of the freshly pre- 
determined .Bfc, by the (m,n) = (—1,0) item of equation (|20p . Due to the symmetry of our 
ansatz for we could have also used (m, n) = (0, ±1) or (m, n) = (0,-1), with the same result 
of course. 



6 Nonequivalent isospectral Hermitian Hamiltonians f) 



In principle, there exist many non-equivalent factorizations of a given quasilocal metric with 
(2R— 1) > 3 diagonals. In our illustrative example we may consider, in general, the superposition 



9 = @(ai,a 2 , ■ ■ ■ ,a R ) 
and factorize G = T 0. 



ai6i + a 2 02 H 1- OR©ii 



(21) 



6.1 R = 1 and the diagonal Dyson mappings 17 

By far the most popular reconstruction of the Dyson operator Q, from a given metric O = Q^Q 
is based on an ad hoc assumption of Hermiticity f2 = & = y/Q (cf., e.g., [22J). For our present 
class of band-matrix models Qr such an assumption does not lead to any problems when the 
metric is diagonal, R = 1. For illustration let us select Oi of equation (I19p . Postulating the 

>t 



diagonality and positivity of the related Dyson map Oi in its factorization ©i 
at the following unique result, 



fijfii we arrive 



v/(l + <?)/(! -<?) 
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knowledge of which enables us to evaluate the related matrix of the Hamiltonian in the physical 
Hilbert space 



2 



2 



(22) 



The off-diagonal elements -y/l - g 2 = -1 + (1 - y/l - g 2 ) = -1 + g 2 /(l + y/T 
composed of the kinetic-energy term (—1) and a positive function of the coupling. 



g 2 ) are 



6.2 R = 2 and a tridiagonal Dyson mapping ft 



We should issue a warning that the usual assumption of Hermiticity of factors fi = fit would 
be rather counterproductive because the Hermitian square root of our sparse matrices © would 
be a non-sparse matrix. The quasilocality property (i.e., a convergence to locality during the 
limiting transition h — ► 0) would be lost for the Dyson map. In opposite direction, when we 
preserve the quasilocality (i.e., the band-matrix structure) of Dyson matrices, we achieve, at 
any nonzero h > 0, a significant simplification of the physical representation of our operators 
of observables in H.^ . 

Without any significant loss of generality the merits of such a strategy may be illustrated 
via the most general tridiagonal metric ©(7) = 2Gi +762- First of all we have to guarantee its 
positivity and invertibility but both these properties are easily shown to be guaranteed inside 
the open interval of 7 G (—1,1). 

Secondly, we have to demonstrate the practical feasibility of the rather difficult extraction of 
a quasilocal factor fi from a given quasilocal = fitfi. For this purpose, let us slightly simplify 
the presentation of the argument and choose 7 — ► — 1 lying on the boundary of the open domain 
of existence of the metric operator (|2ip . In this limiting case our doubly infinite metric acquires 
a rather elementary matrix form, 







2 -1 
-1 2 
-1 



-1 

2 -l-g 
■1-9 2±±§ 
-l-g 



■l-g 
2 -1 
-1 2 -1 

-1 2 



This matrix is most easily factorized into band-matrix Dyson-mapping factors, = fitfi when 
we postulate the following special form of the Dyson-mapping matrix which is asymmetric and 



Cryptohermitian Picture of Scattering Using Quasilocal Metric Operators 



17 



tridiagonal, 





1 -i 








i -i 








i 


-l-g 




n = 




/2<? 2 (l+9) 
V 1-9 








-l-g 


1 

-1 1 

-1 1 



(23) 



During the transition to a measurable coordinate y m appearing as the argument in the physical 
wave function tp^(y m ) S TL^\ the similar quasilocal Dyson mappings enter the scene via finite 
sums 

^ P) (y m ) ■= ^0(y m „x n )^ (F) (x n ). (24) 

n 

In an effective theory as advocated in [25], both the arguments x m and y m may be considered, in 
some sense, observable. Then, equation (|24p introduces just a certain "smearing" of coordinates, 
at the non-vanishing lattice sizes h > at least. The extent of this smearing is proportional to 
the number of diagonals in f2. 

We have to remember that we choose the parameter 7 = 1 which lies, strictly speaking, out 
of the open interval of its admissible values. Due care is needed when working with this option. 
Nevertheless, a part of its undeniable methodical appeal still recurs with the easiness of the 
evaluation of the inverse 



1 1 1 
1 1 
1 


/l-9 2 

V 2g 2 

/l~9 2 
V 2^ 

V w 






I 1~9 




V ^(i+g) 




/1-9 2 

V 2g 2 

/l-9 2 
V 2 9 2 
/l-9 2 

V V 


1 

1 1 
1 1 1 



The extreme 7 = 1 example also offers a very useful guide for transition to the generic \j\ < 1 
where an ansatz for ^(7) may be made with the same sparse-matrix structure. This would lead 
to an efficient factorization recipe based on the use of continued fractions in the definition of 
inverse matrix (at R = 2, cf. its sample in [30]). A generalization of this algorithm to R > 2 
using the so called extended continued fractions also exists (cf. [3T]). 

Let us now return to our R = 2 schematic example with the initial Hamiltonian H of equa- 
tion (|18p and with the tridiagonal metric 6(7) where 7 = 1. Easily we may show, by direct 
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computation, that the related isospectral Hermitian Hamiltonian f) reads 



'•• 2 -1 

-1 2 


-1 




-1 


2-g 2 


-V2<7 2 (l-5 2 ) 


i-g 2 






-V2 9 2 (l-5 2 ) 
1-9 2 


-y/2g*(l-g*) 


-^2 5 2 (l-5 2 ) 

2-g 2 


-l 








-l 


2 -1 
-1 2 '•■ 



We see that it remains sparse but different from its diagonal-metric predecessor (|22|h 



7 Discussion 



Beyond the horizons given by our present sample of "PT-symmetric H we may expect that many 
other models would also profit from the h > approximation treating the coordinate x as 
a discretized quantity Xk = hk. This approach enabled us to construct the (2R — l)-diagonal 
metrics at all R and it also allowed us to postpone the study of the continuous coordinate limit 
h — > to the very end of all the calculations. Let us note now that this limiting transition need 
not necessarily be easy. Sometimes, it may be necessary to replace the simple-minded function 
@(x,x') of two real variables in formula (jlOp by a suitable (e.g., momentum-dependent) operator 
generalization. 

At a fixed level of approximation h > one can usually skip the difficult discussions of the 
continuous coordinate limit. Even then, due care must be paid to the (approximate) discrete 
theory where the values x n of the individual coordinates play the role of arguments in wave 
functions i/)^(x„) G without probabilistic interpretation. The same discretized coordi- 

nates also enter the correspondence-principle-reflecting definitions of the potential and of the 
kinetic-energy operator A of equation (I12j) . All of these operators certainly carry the decisive 
information about the dynamics of the quantum system, provided only that the wave functions 
ip( F \x n ) are mapped, into the physical space 7i.( p \ according to the following diagram, 



ip transferred in TC^^ : 
physics = transparent 
calculations = prohibitively difficult 



simplification 



\ *\ unitary equivalence 



ip calculated in : 
calculations = feasible 
physical meaning = lost 



hermitization 



ip reinterpreted in : 
metric 0^7 = strange 
inner product = unusual 



Physical predictions necessitate a transfer of wave functions from the computation-facilitating 
Hilbert space TL^ to the correct physical Hilbert space 7i^\ The ambiguity of this transfer is 
well illustrated by equation ((Jj) containing the set of free parameters Xj and Kj which reflects 
their presence also in the related Dyson map £1 — r^(Aj, Kj). 

The main consequence of such a departure from the dictum of textbooks is that for a quantized 
point particle our construction and knowledge of wave functions must be complemented by 
a suitable upgrade of Hilbert space, i.e., usually, of the unphysical L 2 (M) := 7i.( F \ In our 
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present models it has to be accompanied by the limiting transition h 
unphysical version of the discrete normalization condition 



^2 ip*(xk)^(xk) 



1 



m 



0. In this way the 



(25) 



k=0 



will be replaced by its continuous limit (cf. equation ([9])) while, similarly, its physical, "stan- 
dardized" counterpart 



^2^2r(x j )Q(x j ,xMx k ) = l in 

j=0 k=0 



(26) 



will be assumed convergent to a double integral formula with h — > (cf. equation (jlOp ). 

In this setting the main message of our present paper can be read as a proposal of transition 
from the diagonal Runge-Kutta metrics (say, of [3 [23]) to their non-diagonal, {2R— l)-diagonal 
sparse-matrix generalizations with R > 2. Under this generalization the latter assumption 
of convergence Q(xj,x' k ) — > Q(x,x') need not always be satisfied since, in general, we have to 
expect the emergence of a more complicated structure in the quasilocal metric = lim Q(xj, x' k ). 

h— >0 

Typically, in a way inspired by the inspection of equation ([7]) we might expect the emergence of 
a complicated dependence of our quasilocal metrics on momenta, etc. 

The study of the criteria distinguishing the metrics with a smooth x-dependence (i.e., with 
operators represented by the mere functions of two variables Q(x,x')) from the other operators 
of metrics without such an elementary representation lies far beyond the scope of our present 
paper. Only a few comments may be added. 

In the first one we repeat that the physical probabilistic interpretation of the system requires 
the introduction of the "standard" Hilbert space ^ L 2 (]R) where the time-evolution becomes 
unitary. Although this space may be defined as spanned by the same functions of x, their norm 
must be defined differently. This postulate degenerates back to the standard textbook scenario 
when one returns to the Dirac's local metric, Q(x, x') — > 5(x — x'). Vice versa, for non-Hermitian 
H 7^ with real spectra one may, sometimes, succeed in constructing the corresponding ad 
hoc metric kernel Q(x,x') ^ 5{x — x'). 

Our second comment will emphasize that our work with difference-equation approximants 
has been motivated by the tedious nature of some alternative perturbation-expansion techniques 
to Schrodinger equations, say, of [221 125] . In the conclusion let us mention a particularly inter- 
esting possibility of a change of perspective. It emerged with the publication of paper [32] where 
the first-quantized Klein-Gordon equation has been considered. In this case, the three-Hilbert 
space formulation of quantum mechanics offers certain interesting new possibilities of the inter- 
pretation of the role of individual spaces. The point is that the initial Hamiltonian H ^ 
acquired the physical meaning (i.e., relativistic covariance) in rather than in H 1 ?) . In our 
present language this would lead to the modified theoretical arrangement of our three Hilbert 
spaces, 



in the third space 
textbook physics, f) = h' 
(ambiguity of f) = QHQ^ 1 ) 



\ equivalence 



in the first space 
kinematics ==>• H ^ W 
(relativistic covariance) 



quasi — her mit izat ion 



in the second space 
H = H* = O^H^Q 
(ambiguous O = ) 



New ways towards old problems could be sought /found in this direction. For example, in the 
light of some recently obtained new results on the first quantization of relativistic particles with 
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spin [33], an extension of these studies to a scattering arrangement (e.g., along the lines indicated 
in our present paper) would be a particularly challenging task. 
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